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ABSTRACT 

The  deqenerate  parabolic  system  (1.1)  in  the  introduction,  serves  as  a 
model  for  heat  conduction  in  a  heteroqeneous  medium  consistinq  of  two  com¬ 
ponents.  The  first  component  is  made  up  of  small  pieces  suspended  in  the 
second  component,  and  the  second  component  underqoes  a  chance  of  phase  at  a 
prescribed  temperature.  This  phenomenon  occurs  in  a  mixture  of  qravel  and  wet 
soil  (for  example  meltinq  of  frozen  soil).  Existence  and  uniqueness  results 
of  weak  solutions  of  the  degenerate  parabolic  problem  are  shown  by  employinq 
monotone  operator  theory.  Local  regularity,  such  as  continuity  and  bound¬ 
edness  of  the  solution  is  studied.  A  discussion  is  provided  about  the  mutual 
interplay  of  the  thermodynamic  temperature  (the  temperature  in  the  first 
component)  and  the  conductive  temperature  (the  temperature  in  the  second 
component). 
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SIGNIFICANCE  AND  EXPLANATION 


Consider  the  melting  of  frozen  soil  (a  mixture  of  gravel  and  moisture )r 
confined  in  a  bounded  region  G.  The  variables  of  the  problem  are  the 
temperature  0  in  the  gravel  and  the  temperature  V  in  the  water.  At  the 
interface  separating  the  liquid  from  the  solid  there  is  a  relationship 
resulting  from  balance  of  energy  (Stefan  type  condition).  The  model  results 
in  a  degenerate  parabolic  system  for  0  and  V>.  The  paper  supplies  the 
existence  and  unigueness  of  0  and  in  a  distribution  space.  More 
interesting  for  the  applications  is  the  pointwise  definition  of  0  and  ^  * 
It  is  shown  that  in  face  0  and  are  continuous.  The  interplay  between 
0  and  V  is  also  investigated.  In  particular,  it  is  shown  that  the  liquid 
phase  l<p  >  0]  is  contained  in  the  region  where  the  rocks  have  positive 
temperature,  and  that  such  a  region  expands  with  time. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  FREE-BOUNDARY  PROBLEM  FOR  A  DEGENERATE  PARABOLIC  SYSTEM 


Emmanuel e  DiBenedetto ^  and  R.  E.  Showalter*^ 

1.*  Introduction.  Wc  shall  be  concerned  with  the  first  Initial -boundary -value  problem 
for  non-negative  solutions  of  a  system  of  nonlinear  partial  differential  equations  of 
the  form 


(1.1) 


a(0(x,  t>>  +h(9(x,  t)-«p(x,t))  3  fjfr.t)  , 

^  P(«K*.t))-A<|>(x,  t)  +h(«|>(x,t)-9(x,t))  3  f2(x,  t) 


and  a  related  free-boundary  problem  of  Stefan  type.  Here  A  denotes  the  Laplacean 
in  the  spatial  variable  x  e  Ra,  h  >  0,  and  the  pair  a ,  3  are  maximal  monotone 
graphs  in  RxR.  If  the  first  equation  were  to  contain  the  term  "-kAC(x,t)M  with 
k  >  0,  then  the  system  (1.1)  would  be  parabolic.  The  situation  we  consider  here 
with  k-0  Is  accordingly  a  degenerate  parabolic  system. 

Although  the  system  (1.1)  with  the  (possibly  multi-valued)  nonlinear  monotone 
graphs  a  ,  3  Is  of  mathematical  Interest  In  Its  own  right,  we  present  In  Section  2 
an  extensive  discussion  of  how  such  a  system  arises  as  a  model  of  heat  conduction  In 
a  composite  material  consisting  of  two  components  in  which  a  change  of  phase  occurs 
In  the  second  component.  This  model  Is  described  by  (1.1)  with  3  obtained  In  the 
special  form  3  (x)  ■  bx  +  LH(x)  where  b  >  0,  L  is  the  latent  heat  of  fusion  and 
H(-)  Is  the  multi-valued  Heaviside  step  function.  Certain  models  of  diffusion 
through  fractured  porous  media  lead  to  the  same  system. 

Our  results  on  (1.1)  are  organised  as  follows.  In  Section  3  we  prove  that  the 
first  initial -bound ary -value  problem  for  (1.1)  Is  well-posed  when  the  data  satisfies 
certain  Integrablllty  conditions  and  3  is  defined  everywhere  on  R.  If  the  data 
Is  non-negative  then  the  solution  Is  likewise  non-negative;  this  property  Is 
essential  for  the  model  problem  discussed  In  Section  2.  We  make  extensive  use  of 
the  theory  of  maximal  monotone  operators  In  Hilbert  space  to  which  we  refer  to 
(2,31. _ 
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Certain  properties  of  the  solution  are  obtained  when  ve  restrict  attention  to 
the  case  where  a  has  a  lower  linear  bound  and  p  *  bl  +LH  as  In  the  model  problem. 
In  Section  4  we  show  that  if  the  data  Is  essentially  bounded  then  the  solution  of 
(1.1)  is  essentially  bounded.  Additional  conditions  on  the  data  are  shown  to  Imply 
that  0  and  <p  are  continuous.  In  order  to  obtain  these  regularity  results  we 
found  It  very  useful  to  treat  the  problem  as  an  equation  of  evolution  rather  than 
to  have  formulated  it  as  a  variational  Inequality. 

In  Section  5  we  exhibit  an  explicit  lover  bound  on  the  first  component  of  the 
solution  of  (1.1).  This  implies  that  the  set  of  points  where  this  function  is 
positive  (the  positivity  set)  is  non -decreasing  with  time.  Finally  we  show  that  the 
positivity  set  of  the  first  component  contains  that  of  the  second  component,  and  an 
example  Is  given  to  show  this  containment  may  be  proper. 

2 .  Diffusion  In  Heterogeneous  Media 

Ve  begin  with  a  mathematical  description  of  diffusion  processes  within  a  medium 
consisting  of  two  components.  A  fundamental  assumption  is  that  the  first  component 
occurs  In  small  Isolated  parts  that  are  suspended  in  the  second  component.  This 
situation  arises  in  thermal  conduction  through  rocky  soil,  since  the  rocks  are  iso¬ 
lated  within  the  soil.  It  also  occurs  In  the  diffusion  of  liquid  or  ©as  through  a 
porous  media  that  has  been  fractured,  since  the  blocks  of  the  medium  are  isolated 
from  one  another  by  the  system  of  fissures.  Next  we  shall  formulate  a  free-boundary 
problem  of  Stefan  type  that  results  from  a  change  of  phase  In  the  second  component 
of  the  medium.  This  arises  In  the  model  of  heat  conduction  through  the  moisture  in 
rocky  soil  since  the  soil  moisture  may  freeze  or  thaw  with  a  corresponding  release  of 
latent  heat;  there  Is  no  moisture  in  the  rocks.  If  ve  consider  diffusion  In  a 
fractured  medium  In  which  the  system  of  fissures  is  only  partially  saturated  then 


>  /  i  ri  r*i , 


wa  can  think  of  the  fissures  as  containing  holes  In  vhich  a  certain  amount  (per 
voltm)  of  the  liquid  or  gas  Is  trapped  and  can  no  longer  take  part  In  the  diffusion. 
Such  a  diffusion  process  is  formally  equivalent  to  the  preceding  heat  conduction  prob¬ 
lem.  Finally,  we  give  a  weak  formulation  of  this  one-phase  Stefan  problem  for  a  two- 
component  medium. 

Consider  the  conduction  of  heat  through  a  heterogeneous  medium  Gcr“  con¬ 
sisting  of  two  components.  As  our  f^rst  model  for  such  a  process  we  take  the  system 
of  equations  In  the  region  Cl  zGx  (0,«) 


(2.1) 


ag^-kAQ+h (9-<p)  ■  f ^ 


b  5*- A*  +  h(<p-0)  “  f2 


where  9  and  9  are  the  temperatures  in  the  first  and  second  components,  respect¬ 
ively.  Each  Is  a  function  of  position  x  e  G  and  time  t  >  0  and  Is  obtained  at 
a  point  x  by  averaging  the  temperature  of  the  corresponding  component  In  a  neigh¬ 
borhood  which  contains  a  sufficiently  large  number  of  pieces  of  both  components. 

The  constants  a  ,  b  are  specific  heats  of  the  respective  components,  k  Is  the 
conductivity  of  the  first,  the  conductivity  of  the  second  component  Is  normalized 
to  unity,  and  the  positive  number  h  is  related  to  the  surface  area  common  to  the 
two  components.  Thus,  h  Is  a  measure  of  the  homogeneity  of  the  material.  The 
system  (2.1)  is  Just  a  pair  of  classical  heat  conduction  equations  together  with  a 
linear  coupling  to  model  the  simplest  exchange  between  components.  Our  basic 
assumption  that  the  first  component  occurs  in  small  parts  Isolated  by  the  second 
component  Implies  that  k*0  in  (2.1).  That  is,  the  particles  of  the  first  com¬ 
ponent  may  store  h*at  (a  >  0)  or  may  exchange  with  the  surrounding  second  com- 
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ponent  (h  >  0),  but  they  cannot  pass  heat  directly  to  other  first -component  par¬ 
ticles  (k*0).  This  Is  the  sense  In  which  the  system  (2.1)  Is  degenerate  parabolic. 

Suppose  there  Is  a  solid-liquid  phase  change  In  the  second  component  at  the 
temperature  <p»0.  We  consider  here  the  (one-phase)  situation  wherein  >  0  every¬ 

where.  The  region  fl  Is  separated  Into  a  conducting  region  fl+  where  <p  >  0  and 
a  non-conducting  region  where  <p*0;  these  correspond  to  completely  melted  and 
partially  frozen  parts,  respectively.  We  need  not  assume  consists  exclusively 

of  ice  but  only  that  It  Is  a  mixture  of  Ice  and  water  in  thermal  equilibrium  at  the 

we  Introduce  the  fraction  of  water, 
is  the  maximal  monotone  Heaviside 
and  H(s)«0  for  s  <  0.  The  two 
S(t).  If  we  let  n  be  the  unit 
the  speed  of  S  (t)  along  n,  then 

and  (1-0  is  the  fraction  of 

ice.  Moreover,  If  H •  (M.  ,  N.  , . . .,  N  ,  N  )  denotes  the  unit  normal  on  the  Interface 

l  z  m  t 

S-U  (  (S  (t), t) ),  we  find  that  (2.2)  Is  equivalent  to 
(2.3)  V  ’  ^l . V  *  LVW)  • 

Each  of  (2.2)  and  (2.3)  is  called  the  Interface  or  free-boundary  condition. 

It  Is  worthwhile  to  recall  the  simple  experiment  In  which  one  applies  a  uni¬ 
form  heat  source  of  Intensity  F  to  a  unit  volume  of  Ice  at  temperature  <p*0. 


melting  temperature.  At  each  point  (x, t)  of  ft 
g(x,  t);  note  that  $  c  H(<p)  In  fl  where  H(-) 
graph  given  by  H(s)  -  1  for  s  >  0,  H(0)  *  [0,1], 
regions  are  separated  at  time  t  by  an  Interface 
normal  on  S(t)  directed  towards  and  V  be 

we  obtain  the  condition 


(2.2) 


-LV(1-|)  on  S(t) 


dm  -*» 

Where  *  Vx<P  ■  n  Is  the  heat  flux  across  S(t) 


The  temperature  remains  at  zero  until  L  units  of  heat  have  been  added.  During 
this  period  there  is  a  fraction  l  of  water  coexisting  with  the  ice  and  |  in¬ 
creases  at  a  constant  rate  F/L.  When  all  the  ice  has  melted,  i  *  1  and  the 
temperature  9  begins  to  rise  at  the  rate  F/b.  The  constants  L  and  b  are 
the  latent  heat  and  specific  heat,  respectively.  We  can  summarize  the  above  by 
stating  that  the  rate  of  Increase  of  the  internal  energy  or  enthalpy  v*b<p  +  L| 

Is  given  by  F.  Later  we  shall  see  that  not  only  is  enthalpy  the  natural  variable 
to  determine  the  state  of  the  process  but  that  it  Is  mathematically  the  proper 
variable  by  which  to  describe  the  evolution  of  the  process. 

Ve  can  now  formulate  our  problem.  With  the  notation  above  we  seek  a  triple 
of  non-negative  real-valued  functions  0  ,  9  ,  |  on  ft  which  satisfy  the  following: 


(2.4) 

•  +  M0-«p)  - 

(2.5) 

and  5  e  H(<p)  , 

In 

ft 

(2.6) 

b  If-A*  +  h(<H>)  -  f2 

In 

(2.7) 

L  §f-  f2  +  hO 

in 

ft0 

(2.8) 

|*  +  LV(1-|)  -  0 

on 

s 

(2.9) 

9*0 

on 

dGx  (0,») 

(2.10) 

0(x,  0)  -  #0(x)  ,  <f>(x,  0)  -  <P0(x) 

and  t(x,0)-{0(x) 

on 

G  . 

The  data  consists  of  the  strictly  positive  numbers  a,b,h,L  and  the  non-negative 
functions  f^  ,  f^  on  ft  and  on  G  for  which  we  assume  ^Q(x)  c  H(^(x) ) 
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for  all  x  c  G.  As  before,  we  have  set  t)  e  ft:  <p(x,  t)  >  0)  and 

ftQ»((x,t)  e  ft;  <p(x,t)*0}.  The  unknown  Interface  S  between  ft+  and  ft^  is 
the  primary  difficulty  in  the  problem. 

It  is  appropriate  to  obtain  a  weak  formulation  of  the  problem  (2.4)- (2. 10) . 
Thia  Is  necessary  even  with  smooth  data  because  the  free  boundary  S  may  vary  in 
a  discontinuous  manner  and  it  is  also  convenient  because  It  casts  the  problem  in¬ 
to  the  form  of  an  evolution  equation  in  Hilbert  space.  Thus  we  first  compute 
in  the  sense  of  distributions  on  ft.  For  each  test  function  *  €  c”<n) 
ve  obtain 


-  -J  c*  |^+«(>w) 

(^•(b9+L)  -  Aip)i  + 

+1  -L 

n0  ^0 

<b  |*-Art*  +  f  a  + 

+  “o 

We  have  assumed  the  interface  S  and  the  restrictions  of  q>  and  |  to  ft+  and 
to  ft^  are  sufficiently  smooth  to  apply  Gauss'  theorem.  This  calculation  shows 
that 

^r- +  h(<p-G)  -  f2 

in  0* (ft)  if  and  only  if  (2.6), (2.7)  and  (2.8)  hold.  From  these  remarks  we  obtain 
the  following  weak  or  generalized  formulation  of  the  two-component  Stefan  problem: 
given  T  >  0  and  the  non-negative  functions  f^  ,  f2  on  ft  and  0Q ,  on  G 

with  4q  €  H(<Pq)*  find  a  non-negative  triple  of  functions  which  satisfy 


t  (V<p-  (Nt . N)  4-LU-DN  _>* 

l  m  t 


T  {(*  V<p-<P  Vi)-(N.,...,N !  ).(b«p+L)N  *} 

•V  1  °  1 
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(2.11) 

0  e  H^O.TjL^C)),  v  e 

L2(0,T;hJ(G)),  v  e  Hl(0,T;H"l(G))  , 

(2.12) 

a  +  h  (9-<p)  »  f  x  In 

L2(0,T;L2(G))  , 

(2.13) 

fj-  dip  +  h  (<p-0)  =  f2  in  L2(0,T;H’l(G))  , 

(2.14) 

<p  -  CbI+LH)_1(v)  in 

L2<0,T;H*(C))  , 

(2.15) 

9(0)  =  0Q  and  v(0)  » 

b<<>0  +  L5q  in  L2(G)  . 

Certainly  a  smooth  solution  of  (2. 11) -(2. 15)  for  which  the  level  set  S  is  a 
smooth  manifold  necessarily  satisfies  (2 .4) -(2. 10) . 

Remarks.  The  condition  b  >  0  arises  later  In  the  discussion  of  properties  of 
solutions  so  we  briefly  indicate  the  significance  of  this  assumption.  The  constant 
b  is  a  measure  of  the  storage  capacity  of  the  second  component  and  it  depends  on 
the  type  of  material  and  also  the  percentage  present  in  the  second  component. 
Similarly  the  constant  L  is  determined  by  the  type  and  percentage  of  this  material 
in  the  second  component  of  the  medium.  The  essential  interest  here  is  in  the  change 
of  phase  phenomenon  so  we  are  concerned  with  the  case  of  a  sufficient  percentage  of 
the  second  component  material  being  present  to  permit  L  >  0.  The  corresponding 
physically  significant  case  is  that  of  b  >  0;  otherwise  we  would  be  considering 
the  unlikely  case  of  a  material  with  positive  latent  heat  of  fusion  but  with  null 
heat  capacity.  Nevertheless,  most  of  our  results  to  follow  are  obtained  from  the 
weaker  assumptions  that  1.  >  0  and  b  >  0. 

The  type  of  the  problem  we  have  called  degenerate  parabolic.  In  the  system 
of  partial  differential  equations  (2.1)  with  k-  0  it  is  of  interest  to  consider 
the  case  of  b*0  (1,7];  one  can  then  reduce  it  to  the  single  partial  differential 
equation 

(a/h)A<p)  -Aq>  =  f^  +  +  (n/h) 
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which  It  of  pseudo-par abol ic  type  (6,22].  This  is  distinctly  not  the  case  for  the 

free-boundary  problem  considered  here.  An  elimination  of  0  from  (2. 12), (2.13) 

leads  to  the  evolution  equation 

.2  df 

(2.16)  (a/h)  — -  +  —  (v  +  a<p- (a/h)Atp)  -Aq>  *  ^  +  f2  +  • 

dt 

The  pair  of  equations  (2 .14),  (2. 16)  gives  an  equation  for  q>  which  is  of  second 
order  in  tine -derivatives,  definitely  not  pseudo-parabolic  unless  both  b*0  and 
L  ■  0.  Thus,  even  in  the  case  of  b  *  0  where  the  local  description  of  the  problem 
contains  a  pseudo -parabolic  equation  (cf.  (2.4)  and  (2.6)  in  ft+),  the  free -boundary 
problem  with  L  >  0  is  not  of  this  type.  Problems  where  the  phase  is  determined 
by  the  first  component  can  be  pseudo-parabolic;  see  [10,18]. 

3.  Existence  and  Uniqueness  of  the  Weak  Solution 

We  shall  prove  that  the  weak  formulation  of  the  Stefan  problem  (2. 11) -(2. 15) 
is  well-posed.  This  will  be  achieved  by  showing  that  the  problem  corresponds  to 
an  evolution  equation  whose  solutions  are  determined  by  a  nonlinear  semigroup  of 
contractions  and  that  the  generator  of  this  semigroup  is  a  subgradient  operacor* 

The  existence  and  uniqueness  of  a  generalized  solution  of  the  Stefan  problem 
la  contained  In  the  following. 

Theorem  1^.  Let  a  and  0  be  maximal  monotone  graphs  on  H  x  K  and  let  J  and 
k  be  proper  convex  lower -semi -continuous  functions  whose  subgradients  are  given 
by  d'J  *  a  ^  and  dk  *  0  Assume 

u0  €  L2(C>,  j(u0>  €  ll  (G),  v0  e  L1  (G)  0  H_1  (G) ,  k(v0>  c  L 1  (G)  , 
e  L2(0,T;L2(G)),  f2  e  L2(0,T;H*l(G))  , 
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and  that  the  domain  of  0  la  equal  to  R.  Then  there  exists  a  unique  quadruple  of 
functions  which  satisfy 

u  e  Hl  (0,T;L2(G)),  v  e  Hl(0,T;H-1(G)) 


(3.1) 

0  e  L2  (0,T;L2  (G)),  <j>  e 

L2(0,T;hJ(G)) 

(3.2) 

5T  +  “<9-9)  -  fx 

in  L2  (0,  T  ;L2  (G)  )  , 

(3.3) 

+  Hv-9)  -  f2 

in  L2(0,T;H*l(G))  , 

(3.4) 

u  c  a(0),  v  «  P(<p) 

a.e.  in  ft  , 

(3.5) 

u(0)  -  uQ  ,  v(0)  -  v0 

a.e.  In  G  . 

2  1 

(a)  If  in  addition  there  is  a  pair  0^  e  L  (G),  9^  e  HQ(G)  for  which  c  cl(9q) 
and  vQ  e  0(9Q)  a.e.  in  G,  and  if  f^  e  R^fO^L2  (G)),  f2  e  (0,T;H~*  (G)),  then 

57  £  L"(0tT;L2(G)),  e  X.“(0,T;H-l(G)),  0  e  L*(0,T;Z.2(G)),  <p  e  l”(0,T;hJ(G))  . 

(b)  If  In  addition  a(0)  3  0,  0(0)  3  0  ami  each  of  the  functions  f^  ,  »  u^  and 

Is  non-negative,  then  each  of  u,v,  0  and  9  is  non -negative . 

2  1 

Proof:  Let  V  be  the  product  space  L  (G) x  HQ(G)  which  has  the  dual 
V*-L2(G)  xH*X(G).  Define  B  e  st(V,V*)  by 

Bu(v)-J  {hCu^-UjHVj-vp  +7u2-Vv2)  ,  u-  [t^  ,  u2]  ,  v-  [v2  ,  Vj]  e  V  . 

G 

1/2  * 

Renorm  V  with  the  equivalent  norm  (Bu(u>)  so  that  B:  V  *►  V  is  the  corres¬ 
ponding  Riesz  isomorphism  of  the  Hilbert  space  V  onto  its  dual.  Note  that  B  is 
given  in  0* (G)  In  the  form 
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Bduj  ,  u2I)  -  [h(u1-u2>  ,  h (u2*uj)  -Au2]  ,  ,  u2J  €  V 


* 

Next  we  consider  the  function  J:  V  *  R  U  (  +  »}  defined  by 

f  J  OfojJ+ktUj))  if  ui  e  l2»  Kuj)  e  Ll  ,  u2  e  L^H  1  ,  k(u2>  e  L1  , 


J(u) 


+  »  otherwise  , 


“  *  [Uj  ,  u2]  e  V 


From  [3,  pp.  115,123)  we  find  that  J  is  a  proper,  convex  and  lower-semi- continuous 
function  on  V*.  Furthermore  the  subgradient  of  J  is  determined  as  follows: 

g  e  dJ(u)  with  g  «  C »  g2l  and  «  ■  ,  u2)  in  V* 

if  and  only  if  for  some  v=  [v^  ,  v^]  e  V  we  have 

g  *  B(v)  and  v^  e  dJC^),  v^  c  dk(u2)  a.e.  in  G  . 

These  computations  are  imnediate  from  the  corresponding  results  of  [3]  on  the  com- 
* 

ponents  of  V  . 

It  is  useful  to  characterise  dJ  explicitly  as  a  composition  of  operators  in 
£)' (G).  Thus  we  define  A;  V  >V  by  /»=fa,£j:  that  is, 

* 

u  €  A(v)  with  u  *  [u^  ,  u^]  e  V  and  v  *=  [v^  ,  v2J  e  V 

if  and  only  if  u^  e  a(v^)  and  u2  €  3(v2)  a.e,  in  G.  (Note  that  A  *  is  the 

★  irk 

subgradient  of  J  computed  from  the  Banach  space  V  to  its  dual  V  *V  and  A 

is  the  corresponding  subgradient  of  the  conjugate  of  J  [11).)  From  the  computations 

above  we  have  the  representation  dJ  =  B  •  A  *  as  desired. 
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Given  the  subgradient  operator  cU  on  the  Hilbert  space  V*,  it  is  well-known 
[2,3]  that  the  initial-value  problem 

+  dJ(v(t))  )f(t)  ,  a.e.  te  (0,T] 

(3.6) 

w(°)  =  w0 

^  it  _  2  * 

has  a  unique  solution  we  H  (0,T;V  )  whenever  w^  e  dotn(J)  and  f  e  L  (0,  T;V  ) 
are  given.  Furthermore,  if  w^  e  dom(dJ)  and  f  e  r1(0,T;V*)  then  this  solution 

dW  DO  ★ 

satisfies  e  L  (0,T;V  ).  These  remarks,  with  the  identifications 
w(t)  -  (u(t),v(t)]  e  V*,  W0=[U0.V0J>  f(t)  *  [f1(t),f2(t)l  and 

[9(t),<p(t)]  =  (f  (t)-w'  (t) )  €  (w(t))  , 

show  that  (3.6)  is  equivalent  to  (3.1) -(3.5)  and  thereby  establish  all  but  (b)  of 
Theorem  1.  For  the  proof  of  (b)  we  first  change  the  data  as  follows:  (i)  Set 
J(s) *  J(0)  for  s  <  0  and  leave  the  values  as  originally  given  for  s  >  0;  thus 

dom(a)  c  [0,  +»).  (ii)  Add  to  3(s)  the  quantity  s  for  those  s  <  0  and  leave 

the  values  as  originally  given  for  s  >  0;  thus  3  is  strictly  monotone  on  C-®,0]. 

Since  Uq  is  non-negative  the  hypotheses  of  Theorem  1  still  hold  so  there  is  exactly 

one  solution  of  (3.1)- (3.5)  with  the  modified  data;  we  denote  it  by  u,v,9,<p  as  be¬ 
fore.  Since  the  domain  of  a  contains  only  non-negative  numbers,  it  follows  that 

9  >  0.  Our  plan  is  to  show  the  remaining  three  functions  are  non-negative. 

Next  ve  consider  the  equation  (3.3)  written  with  right  side  hQ  + f 2  and  initial 
condition  v^  being  non-negative.  This  equation  is  of  independent  interest. 

Lensna  1.  Let  A  =  h  -  A  be  the  indicated  Riesz  map  of  the  Hilbert  space  Hq(G)  onto 

its  dual,  H  *(G),  and  let  H  *(G)  have  the  scalar-product  corresponding  to  A. 


12. 


Let  y  be  a  maximal  monotone  graph  on  R  x  R  which  contains  the  origin  and  whose 
range  la  all  of  R. 

(a)  The  operator  A°y  is  maximal  monotone  on  H  ^(G)  with  the  domain 

{v  c  H  l0L*:  there  is  a  <p  e  H^(G)  with  ip(x)  e  y(v(x))  a.e,  x  e  G} . 

(b)  If  cs  (f  c  H_1(x):  f  >  0},  then  [I  +M  •  y]"l(C)  c  C  for  every  *  >  0. 

Proof:  Part  (a)  follows  from  Theorem  17  of  (3)  where  it  is  shown  that  A  •  y  is  a 

•ubgradlent  on  To  verify  (b),  let  (I  +M  0  y)  (v)  “  f  in  H*1(G)  with 

f  >  0.  By  truncation  and  regularization  we  obtain  a  sequence  f^  e  L1  (G)  fl  C°°(G)  0  C 

with  f^  *  f  in  H  *(G).  Since  A  <»  Y  is  maximal  monotone,  the  corresponding 

sequence  v^hJi+Moy)  *f  converges  to  v  in  H  *(G).  From  Proposition  5  of  [5] 

It  follows  that  each  v  e  C.  so  we  have  v  e  C. 
n 

Let  F  be  a  maximal  monotone  operator  on  a  Hilbert  space  H;  let  v^  e  dom(F) 
and  f  c  L*(0,T;H).  Then  there  exists  a  unique  weak  solution  of  the  initial-value 
problem  (2,  p.  64] 

(3.7)  jj£  +  F(v)  3  f  on  [0,T]  ,  v(0)  -  vQ  . 

By  a  weak  solution  ve  mean  a  uniform  limit  of  strong  solutions  v^  corresponding 
to  data  v®  and  with  v“  *  vQ  and  f  *  f  in  H  and  L1(0,T;H),  respectively. 

This  existence  result  Is  proved  by  choosing  the  sequences  above  with  each  v”  c  dom(F) 
and  each  f  a  step-function  with  values  from  the  range  of  f,  [2,  p.  65]. 

Lemma  2.  Let  C  be  a  closed  cone  in  H.  If  vQ  e  Cf  f(t)  e  C  for  all  t  e  [0,T], 
and  if  [I+%F]  *(C)  c  c  for  all  >»  >  0,  then  the  weak  solution  v  of  (3.7) 
satisfies  v(t)  e  C  for  all  t  e  (0,T]. 


Proof:  By  the  preceding  remarks  It  suffices  to  consider  the  case  of  e  dom(F) 
and  a  step-function  f  given  on  a  partition  0  «  a^  <  a^  a^  *  T  by  f  »  e  C 

on  IaA  j  »  *j)*  The  solution  is  given  inductively  by  v(0)  -v^  and 
v(t)  -S1(t-a1)v(at)  on  fa^  ,  a^]  where  is  the  semigroup  generated  by 

-(F-yj).  By  [2,  Proposition  4.5]  it  suffices  to  show  [I+A(F-y^)]  l(C)  c  C,  for 
then  we  have  v(t)  e  C  for  all  t  e  (0,T].  Thus,  let  x*  (I+MF-y^)]  \  with 
y  €  C  and  A  >  0.  It  follows  directly  that  x«  (I  +>vF)~*  Qvy^  +  y) .  Since 
>yt+y  c  C  we  have  x  e  C  and  we  are  done. 

To  obtain  v  >  0  In  (3.3)  we  apply  Lemma  1  with  y  *  0  *  and  then  apply  Lemna  2 
with  F*A*y.  Since  3  is  strictly  monotone  on  (-«,0]  it  follows  from  (3.4) 
that  9  >  0.  Finally,  writing  u  as  the  sum  of  its  positive  and  negative  parts, 
u -  u+  - u”,  we  obtain  from  (3.2) 


(1/2)  4~  f  (u-)2  -  ,  -  h(u",«)  2  -  (u*f h<p  +  f.)  2 

dtJG  dt  L2(G)  L2(G)  1L2(G) 

Since  u  and  9  have  the  same  sign  and  tup  +  f^  is  non-negative,  the  right  side 
la  non-positive  so  u”«0.  Thus  all  four  of  u,v,  9,9  are  non-negatlve.  It  follows 
that  this  quadruple  is  a  solution  of  the  original  problem  without  the  modified  data. 
By  uniqueness  this  is  the  solution  of  the  original  problem  and  (b)  1$  established. 

Remarks .  The  essential  point  in  the  first  part  of  the  proof  of  Theorem  1  is  to  re¬ 
duce  the  problem  to  the  evolution  equation  (3.6)  whose  solution  is  the  pair 
(u(t)yV(t)]  of  "enthalpy"  functions  associated  with  the  weak  solution.  It  is 

this  sense  in  which  enthalpy  Is  the  natural  variable  for  the  problem. 

2  2 

Por  the  special  case  of  f^  e  L  (0, T;L  (G))  we  can  give  an  alternate  proof 

of  part  (b)  of  Theorem  1  as  follows.  Approximate  p  by  a  smooth  3  for  which  the 
r  n 


9 

corresponding  solutions  [uq  ,  v^]  can  be  shown  to  be  non-negative  by  direct  ^-esti¬ 
mates  on  (3,2)  and  (3.3).  Then  using  methods  of  [9]  we  can  let  n  to  obtain  the 

non-negatlvlty  of  (u,v].  However  the  proof  given  above  permits  the  more  general 
data  of  the  existence  result,  and  we  also  obtain  the  corresponding  well-known  non- 
negatlvlty  result  for  the  abstract  porous  media  equation 

JT  +  (h^)Y(v)  3  f(t) 

in  H  *(G),  where  h  >  0  and  y  Is  maximal  monotone.  We  could  not  find  this  re¬ 
sult  in  the  literature. 

The  evolution  equation  (3.6)  Is  of  the  form 

+  A_l(w(t»  »  f(t) 

where  is  positive  self-adjoint  and  A*1  is  maximal  monotone  from  a  Hilbert 

Space  to  its  dual.  Various  generalisations  and  related  equations  have  been  discussed 
in  [4, 6,  9, 10, 15, 16, 20, 21] . 

4 .  Boundedness  and  Continuity  of  the  Weak  Solution 

We  shall  prove  that  the  ”  temperatures”  9  and  in  the  weak  solution  are 
bounded  when  the  data  in  the  problem  is  bounded.  We  also  give  sufficient  conditions 
for  9  and  <p  to  be  continuous.  These  results  are  obtained  In  the  following 
special  case  of  Theorem  1  which  contains  the  weak  formulation  of  the  one-phase  two- 
component  Stefan  problem. 


Theorem  2.  In  addition  to  the  conditions  of  Theorem  1(b)  we  assume  the  following: 
(1)  there  is  a  number  a  >  0  such  that  r  >  as  for  all  r  e  a(s) ;  (ii)  the  maxi¬ 
mal  monotone  p  is  of  the  form  p-bl  +  LH  where  I  is  the  identity,  H  is  the 
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Heaviside  graph  and  both  b  and  L  are  non-negative;  (iii)  the  initial  data  and 
forcing  terms  are  essentially  bounded:  ,  v^  €  L°°(G)  and  f,  ,  e  t° (ft) .  Then 

the  functions  u,v,0,9  are  bounded  on  ft. 

(a)  If  in  addition  we  have  jr^-r^l  >  als^-s^l  for  all  r^  €  a(a^)  and 

e  and  if  the  functions  and  J^f^(*,t)dt  are  uniformly  (Holder) 

continuous  on  G,  then  0  is  uniformly  (respectively.  Holder)  continuous  on  ft, 

(b)  If  b  >  0  and  is  uniformly  continuous  on  G,  then  9  is  uniformly  con¬ 

tinuous  on  ft . 


Proof:  Let  u,v,0,9  be  the  solution  of  (3.1)-(3.5);  by  assumption  (ii)  we  may 
write  v*b<p  +  L£,  |  €  H(<p),  in  ft.  For  each  t  >  0  we  consider  the  Steklov 
averages 


1  o  *  i  * 

<Pe<t)  «  ~  J  M>(s)ds  ,  0c(t)  =  ±  j  0(s)ds 


t 

t-c 


where  9  and  0  are  extended  as  9(0),  0(0),  respectively,  on  (-e,0).  It  is 

known  that  lim<p  *=9  (e.t.c.)  in  (0, T;H* (G))  [14,  p.85).  By  integrating  (3-3) 

e*0e  0 

over  [t-e,t]  we  obtain 


b9^(t)  +  (L/e)(£  (t)-£  (t-c))  +h9c(t)  -  &9c<t)  =*  h0£(t)  +  (1/c 


f2(s)ds  . 


V/e  shall  apply  this  to  (9(x)-k)+,  integrate  over  (0,  t)  x  G  where  the  superscript- 
plus  denotes  the  positive  part  of  the  indicated  function  in  Hq(G)  and  the  number 
k  is  chosen  by 


k  =max(  ||<(>0||  ;|||u 

0  L  (G)  8  C 


L  (G) 


■bK\ 


+  2|!f„ 


L  (ft) 


I*  (ft) 


))  . 
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and  take  the  limit  as  c  l  0.  To  this  end  we  obtain  the  following: 


ii*  r  r  <p’(x)(<p(z)-k)+  -  (j)iK<p(t).k)+ii22  -(7)iKvk)+ii22 

c>0J0JC  *  2  L  (G)  2  0  L2 


(0) 


and  the  lest  tens  vanishes  since  k  >  ||<pn||  ,  ; 

L  (G) 


(((t)-((t-e))(<p(t)-k)+  >  (|(t)-l)(<i>(t)-k)+  -  0 


since  k  >  0  and  |  (t)  c  H(<p(t>);  and 


It*  J  J  C*c-^c) (f-k)* -  J  J  (h|  (<|>-k)+|2 +hk(cp-k)++ |v(<p-k)+|2) 


e >0  w0  G 


'0  G 


Thus  we  obtain 


(|)ll(«P(t)-k)Y2  +  **(**  f  <9-k)++h  J^f  IC'HoY+f'j  |v<«H0+|2 

2  l/(G)  °0°G  *  n  n  * 


J0  UG 


' 0  G 


<  i.  r'j  •<♦*>*♦£. X  ,(eJT  V(<p(T)’k)+dTdx  * 


'0  G 


0  G  r-c 


This  leads  lamed  lately  to  the  estimate 


<4-i)  k J«  X,(,H0+  • 


0  G 


0  "G 


L  (H)  ~G 


Meat  we  estimate  the  first  term  on  the  right  side  of  (4.1).  Integrate  (3.2) 
over  (0,t)  and  use  (1)  and  0  >0  to  obtain 


(4,2) 


*  t  t  t 

a0(t)  <  u(t)  +h  J  0  -  h  J  <p+Jo  +  uQ 
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From  (4.2)  easily  follows 


•9(t>  <  h  r  GHO+  +  t(hk  +  Ilf  II  )  +  flu  II 

J0  1  L  (fl)  0  L  (G) 


Nov  apply  this  to  (<p-k)  end  integrate  to  obtain 


J  J  Q(«MO+  <  h  J  J  {  (<p(r)-k)  +  J  (<HO+]dxdr 


'0  C 


'0  VG 


+  lt(Kk+||f  ||  w  )  +  (|u  ||  .  ]  J  J  <«H0* 

1  L  (p)  0  L  (G)  J0  JG 


Mote  that 


J  *  J  (  (<P<T)-k)+  J,T(9-k)+)d*dr  <  <|)  (J (<p-k)+||  2  )2  <  (|)  J *||  (v-k)+ 


' 0  G 


so  w  have  from  above 


L  <G) 


L  <G) 


(4.3) 


V+..2 


0  G 


11  2 
l/(G) 


b  (A) 


L  (G)  m0vG 


If  ve  use  (4.3)  In  (4.1)  ve  obtain 


o-U>  <- 1  <v- 1>k  *5nfiii,.  ’X11,-  *s|,!|i,.1  S*S„<**>* 


L  v0  G 


Thus,  if  0  <  t  <  a/2h  then  by  our  choice  of  k  the  right  side  Is  non-posltlve 
and  the  left  side  Is  necessarily  zero. 


A. 
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In  stannary,  we  have  shown  that  with  k  ss  given  above  we  have 


Itooll . 

<  k  , 

L  (G) 

►<*>l  „ 

L  (G) 

«•<*># . 

<  (k  /«)e1/2  , 

L  (C) 

4 

for  a.e.  t  c  [0,a/2h].  The  first  is  Immediate  from  our  preceding  calculations, 
the  second  follows  from  (4.2),  and  the  third  Is  obtained  fro*  (4i?N  md  Grouwall's 
Inequality.  From  the  dependence  of  k  on  the  data  It  is  dee'  -f  h*  Jie  estimates 
(4*4)  on  Gx  (0,  a/2h)  can  be  extended  to  give  a  bound  o n  i**#  9  on  all  of 

fl  In  a  finite  number  of  steps. 

In  order  to  prove  (a)  we  first  consider  the  function* 


it  t  t 

•  (x,t)  -  V  q>(x,»)d»  ,  a  (x,t)  -  y  0(x,svd»  . 
J0  J0 


Integrate  (3.3)  to  see  that  $  is  a  weak  solution  of 

M  r>Z 

bH-^  +  W-h®  +  J  f2  +  L|0  -  14 

k  $(x,0)  *  0  ,  e  H^(G)  ,  0  <  t  <  T 


If  h>  0  then  from  [14,  Theorem  1.1,  p.  419]  we  conclude  that  ft  is  uniformly 
Holder  continuous  on  ft.  If  b  *  0  then  from  [13,  Theorem  14.1,  p.  201]  we  con- 
elude  that  $  is  uniformly  Holder  continuous  on  G,  uniformly  In  t  c  (0,T). 
Since  9  Is  bounded,  <t>  is  trivially  Lipschltz  In  t.  Thus  It  follows  that 
6  la  uniformly  Holder  continuous  on  11. 


i* 

» 


By  Gronwall's  Inequality  It  follows  that  0  has  the  same  modulus  of  continuity, 
op  in  x.  From  (3.2)  follows  the  uniform  Lipschitz  continuity  In  t  of  u  and 
then  the  assumption  in  (a)  shows  that  0  is  uniformly  Lipschitz  in  t .  This 
finishes  the  proof  of  (a). 

The  proof  of  (b)  Is  an  immediate  corollary  of  [8,  Theorem  5.3,  p.  69].  The 
point  is  that  (p  Is  an  essentially  bounded  weak  solution  of 


b  It  +  L  tt-A*  “  +  f: 


i  e  H(<|>)  , 


2  * 

with  c  (H ) .  This  last  Inclusion  follows  from  b  >  0;  see  [9]  or  [14,  p.SOl]. 

Remarks .  The  boundary  bG  of  the  region  G  Is  assumed  to  satisfy  "Condition  A" 
of  (13]  in  both  (a)  and  (b)  of  Theorem  2.  That  is,  there  is  a  pair  of  positive 
masbers  a^  and  such  that  for  any  sphere  with  center  on  dG  of  radius 

A 

r  <  a^  and  for  any  component  of  the  intersection  G^-B^OG  it  follows  that 

mes(Gr)  <  (l-O^)mes  (B^)  .  Without  such  a  restriction  on  the  smoothness  of  the  boundary 
of  G  we  obtain  local  or  interior  continuity  results  as  above. 

It  la  not  known  whether  9  is  continuous  in  ft  in  the  case  b  *  0.  Fortunately 
this  situation  Is  of  no  apparent  interest  for  the  applications. 

Additional  Properties  of  the  Weak  Solution 

Under  rather  general  conditions  on  the  data  in  the  weak  formulation  of  our  prob¬ 
lem  It  follows  that  the  positivity  set  of  the  enthalpy  u  is  increasing  with  time. 
This  is  equal  to  the  positivity  set  of  the  temperature  0  and  we  show  that  it  con¬ 
tains  the  positivity  set  of  the  temperature  An  example  shows  this  containment 

may  be  proper. 

The  preceding  properties  of  the  weak  solution  will  be  obtained  In  part  from  the 
following  comparison  result. 

Lama  3.  Suppose  t^  <  T  and  for  each  t  e  ( t^  ,  T]  we  are  given  a  pair  y^Ct,  •)* 

>2  (tf  •)  of  graphs  on  R  x  R  such  that  y^t,  •)  is  monotone  and 

(for  each  s^  c  y^(t,r)  there  exists  a 
*2  €  Y2(t,r)  for  which  *2  <  \  * 


(5.1) 
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Let  the  pair  of  absolutely  continuous  functions  u^  ,  u2:  [t^.T]  *  R  satisfy 
uj  ^0*  £  “2  *“* 

u|(t)  +Y1(t#u1(t))  a  0,  u^(t)  +v2(t,u2(t))  3  0  , 

for  a.e.  t  c  (t^f  Tj.  Then  u^(t)  <  u2(t)  for  tQ  <  t  <  T. 

Proof:  Suppose  there  Is  a  t2  e  (t^  ,  T]  such  that  u^(t2>  >  u2(t2>.  Define 

tj^lubCt  e  ( t^  ,  T] :  Uj(t)  <  u2(t)}  and  note  that  (t^)  =  u2  (t^)  and 

u^(t)  >  u2(t)  for  all  t  e  Jt^  ,  t^}.  For  each  t  €  [t^  ,  t2]  for  which 
-uj(t)  e  Y1(t,ul(t))  there  is  by  (5.1)  an  s(t)  e  Y2(tt^1(t>)  with  s(t)  < 

For  such  a  t  we  obtain 

(uj(t)-u^(t))  (u^(t)-u2(t))  <  (-s (t) -u *  (t ) ) (t) -u2 (t ) )  | 

and  this  last  quantity  is  non-positive  because  y2 (t, ■)  is  monotone.  Thus  the 
2 

function  (u^ (t)-u2(t))  has  a  non-positive  derivative  on  [t^  ,  t2]»  it  there¬ 
fore  vanishes  on  [t^,  t2J  and  this  contradicts  the  choice  of  t2  . 

Suppose  we  are  in  the  situation  of  Theorem  1(b).  For  a.e.  x  e  G  the  function 

u2(t>£u(x, t)  is  an  absolutely  continuous  solution  of  (3.2).  In  order  to  apply 
Lews  3  to  (3.2),  (3.4),  define 

Y^tjU)  -  (h/a)u+-fx(x,t)  ,  Y2<t»u)  *  ha"1  (u)  -  hq>(x,  t)  -  f^x,  t)  . 

We  will  as s use  a  >  0  and  that 

for  each  r  >  0  there  exists  an  s  e  a~*  (r) 
such  that  as  <  r  . 


(5.2) 
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This  Implies  that  (5.1)  holds  for  our  choice  of  .  Let  he  the  solution  of 


u|(t)  +  (h/a)u1(t)  -  f1(x,t)  ,  0  <  t  <  T 


with  uj(1q)  *u(x,  t^).  From  Lemma  3  we  obtain  the  first  part  of  the  following. 

Theorem  3.  In  addition  to  the  conditions  of  Theorem  1(b)  we  assume  there  Is  an  a  >  0 
for  which  (5.2)  holds.  Then  the  first  component  of  the  solution  of  (3.1)* (3.5)  satisfies 

-7<t-t0)  rt 

(5.3)  u(x,t)  >  e  -u(x,t0)  +  J  e  f^x,  s)da  ,  0  <  tQ  <  t  <  T 

*0 

for  almost  every  x  c  G.  Thus  the  set  S*(u)  s  (x  e  G:  u(x,t)  >  0)  is  Increasing 
with  t.  Furthermore  the  set  S+(u)  =  {(x,t)  e  fl:  u(x,  t)  >  0)  contains  the  Interior 
of  S+(fx)  and 

(5.4)  U(S^(9)i  0  <  t  <  tj)  c  $*  (u)  ,  0  <  tj  <  T  . 

Proof;  The  Inequality  (5.3)  follows  from  the  preceding  remarks  and  it  Immediately 

4-  + 

Implies  the  mono tonicity  of  St(u)  and  the  Inclusion  of  the  interior  of  S  (f^) 

In  S*(u).  We  verify  (5.4).  Let  4  S*  (u),  that  is,  u(x^  ,  t^)  *  0,  so  by 

(5.3)  we  have  u(x^  ,  t)  «  0  for  all  0  <  t  <  t^  .  Thus  du(x^  ,  t)/dt«  0  for 

0  <  t  <  t^  .  From  (5.2)  we  obtain  (see  below)  0(x^  ,  t)  *  0  for  0  <  t  <  t^  so 

(3.2)  implies  ^(x^  ,  t)  «  0  for  0  <  t  <  t^  .  That  is,  x^  4  S*(<p)  for  all 

0  <  t  <  tx . 

Corollary .  In  the  situation  of  Theorem  2(b)  we  have  S+(<p)  c  S+(u) . 


Proof:  Since  <p  is  continuous  this  follows  from  (5.4). 
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Remarks .  The  condition  (5.2)  is  actually  equivalent  to  the  assumption  (i)  in 


Theorem  2:  i  >  as  for  all  (s,r)  e  a.  To  see  this,  note  that  if  (s^  ,  r^)  c  a 
with  rQ  <  «Sq  ,  then  we  can  choose  r^  =  (1/2) (rQ  +  asQ)  >  0  and  from  (5.2)  a 

with  ( ,  r^)  e  a  and  as^  <  .  But  then  <  (1/2)  (rQ/a  +  sQ)  <  and 

Tj  >  ,  contradicting  the  monotonicity  of  a.  Thus,  (5.2)  implies  that  all  of 

a  lies  above  the  graph  of  r=  as. 

As  a  consequence  of  the  above  remark  it  follows  from  (3.4)  that  u(x,t)  >  aO(x,t), 

hence  S+(9)  c  S+(u).  If  in  addition  we  have  a(0)*{0),  then  S+(0)**S+(u). 

In  the  case  of  our  original  problem,  (2 . 11) - (2 . 15),  we  have  S+(<p)  c  S+(0) : 

Thus  9  >  0  in  the  region  where  the  water  is  completely  melted.  The  follow¬ 

ing  example  shows  that  we  do  not  necessarily  have  S+(tp)  =  S+(9)  . 

Example.  Define  9(t)  =  e  t,  q>  »  0  and  v(t)  *  1  -  e  for  t  >  0.  This  triple  of 

non-negative  functions  is  the  solution  of  (2 . 11)- (2 . 15)  with  L  =  a =  h =  1  and 
arbitrary  b  >  0,  f^  *  “  0,  and  9q=1,  =  This  solution  is  independent 

of  x  e  G.  In  the  thermal  conduction  model  of  Section  2,  this  example  corresponds 
to  the  situation  wherein  a  small  amount  of  heat  uniformly  distributed  in  the  first 
component  is  all  absorbed  as  latent  heat  to  convert  the  second  component  from  solid 
ice  to  water  at  temperature  equal  to  zero. 
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